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Using quasiparticle self-consistent GW calculations, we re-examined the electronic structure of 
Sr 2 Ru 04 and SrRuOa. Our calculations show that the correlation effects beyond the conventional 
LDA (local density approximation) and GGA (generalized gradient approximation) are reasonably 
well captured by QSGIU self-energy without any ad hoc parameter or any ambiguity related to 
the double-counting and the downfolding issues. While the spectral weight transfer to the lower 
and upper Hubbard band is not observed, the noticeable bandwidth reduction and effective mass 
enhancement are obtained. Important features in the electronic structures that have been debated 
over the last decades such as the photoemission spectra at around —3 eV in Sr 2 Ru 04 and the half- 
metallicity for SrRuOs are discussed in the light of our QSGIU results and in comparison with the 
previous studies. The promising aspects of QSGW are highlighted as the first-principles calculation 
method to describe the moderately correlated 4d transition metal oxides along with its limitations. 

PACS numbers: 71.20.Be,71.15.-m, 71.18.-|-y 


I. INTRODUCTION 

Since the seminal discovery of unconventional super¬ 
conductivity at < 1 K, Sr 2 Ru 04 (SR0214) has been stud¬ 
ied extensively [T]. The crystal structure of SR0214 is 
of K 2 NiF 4 -type at low temperature as shown in Fig. 
and its normal phase is a paramagnetic metal. Its in¬ 
triguing electronic behavior and superconductivity are 
still a subject of active study [2HS]- SrRuOs (SR0113) is 
a ferromagnetic metal with a transition temperature of 
Tc ~ 160 K. The observed magnetic moment is /!= 1.1- 
1.7 /le/f-u. [TOUM] and the stable structure at low tem¬ 
perature is an orthorhombic perovskite with a GdFeOs- 
type distortion. The electronic structure of SR0113 is lo¬ 
cated near to the half-metallicity as reported by density 
functional theory (DFT) -|- U calculations [UHlT]- The 
thin film SR0113 is of particular interest as a widely- 
used bottom electrode [TBfED] and is intensively studied 
for the possibility of a new field-effect device [2TVI23] . 

Considering their fundamental and technological im¬ 
portance, it is essential to understand correctly the elec¬ 
tronic structures of SR0214 and SR0113. For SR0214, 
the early electronic structure calculations showed 

that the states near the Fermi level (Ep) consist of anti¬ 
bonding combinations of Ru-t 2 g and 0-2p, and the Fermi 
surface is composed of the F-centered electron-like sheets 
(called as /3 and 7 ) and the X-point-centered hole-like 
sheet (called as a) PIHIS]- The LDA (local density ap¬ 
proximation) Fermi surface topology was in good agree¬ 
ment with that of de Haas-van Alphen (dHvA) experi¬ 
ments [301 El] ) but does not seem to be with the angle- 
resolved photoemission spectroscopy (ARPES) |3SJ|50] in 
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the sense that one electron-like (/?) and two hole-like (a 
and 7 ) Fermi surface were observed by ARPES [33130] . 
It however turns out to be a surface effect [37030] later. 
Therefore the overall features of the electronic structure 
can be regarded as being described well by the conven¬ 
tional electronic structure calculation techniques such as 
LDA and GGA (generalized gradient approximation). 

However the details of the electronic structure are still 
not clearly understood. A series of recent studies that 
identify SR0214 as a ‘Hund’s metal’ [2JS| can be one ex¬ 
ample showing that the previous understanding based on 
the conventional electronic structure calculations was not 
enough. Also, there was a debate regarding the detailed 
feature of electronic levels at ~ — 3 eV revealed by X-ray 
photoemission spectroscopy (XPS) [iDM3113711^ . Ac¬ 
cording to XPS and resonant XPS, the indication of lower 
Hubbard band (LHB) of Ru-4d states is observed while 
the LDA calculation predicts that those features should 
be attributed mainly to the oxygen states [37[. The 
experimental effective mass (to*/tolda — 2-5) is also 
significantly larger than the LDA value [301133] . While 
DMET can be a reasonable alternative in this situation 
[H I33M3] , its parameter dependency has caused the dis¬ 
crepancy even among the DMFT results. Eor example, 
the calculations by Pchelkina et al. [45] and by Liebsch 
et al. |46j seem to give different answers to the nature 
of the states at ~ — 3 eV most likely due to their differ¬ 
ent choices of U and J values. Beside the ambiguity in 
the double-counting term in DFT-t-DMFT, this param¬ 
eter dependency is a well known problem for the first- 
principles description of the correlated electron systems. 
It is also noted that another well-established technique, 
DET-I-H, is not useful for paramagnetic SR0214 because 
of its Hartree-Eock nature which always prefers the mag¬ 
netic solution [33] . 

Eor SR0113 the situation is similar with SR0214; 
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while the overall features can be described by the conven¬ 
tional methods, the details are not clearly understood. 
Note that DFT-l-17 can be used for this case as the 
ground state of SR0113 is (ferro) magnetic. For the bulk 
phase, LDA, LDA-I-C/, and self-interaction correction 
(SIC) give the correct ferromagnetic solution although 
the calculated moment shows some deviations [m [SDi- 
15^ . For the thin film SR0113, however, the experimen¬ 
tally observed metal-to-insulator transition (MIT) and 
the ferromagnetic-to-nonmagnetic transition as a func¬ 
tion of layer thickness are not consistently reproduced 
by these techniques [MEM?]. For example, the single¬ 
layer SROI13 is predicted to be either ferromagnetic or 
nonmagnetic depending on the choice of the exchange- 
correlation functional and the U values |5?1 [55H57] . As 
in the case of SR02I4, DMFT has made important con¬ 
tributions for SROII3 [21 m |58l [59], but it is still not 
completely satisfactory because of the limitations such 
as the parameter dependency. Also, the details of the 
spectroscopic data of SR0113 needs further clarifications 

[Ml]. 

In this paper, we re-investigate the electronic struc¬ 
ture of these two classical Ad transition-metal oxide sys¬ 
tems by using quasiparticle self-consistent GW (QSGkF) 
[55H57] method. To the best of our knowledge, this 
well-established calculation technique has never been ap¬ 
plied to the ruthenates while GW self-energy can be ex¬ 
pected to give a good description for the metallic and 
weakly correlated systems. Despite of its limitation to 
take the Hubbard-like on-site Coulomb interactions into 
account, QSGVF has a distinctive advantage as a fully 
self-consistent ‘parameter-free’ technique. Due to the 
recent experimental progress in making ‘complex ox¬ 
ide’ structures such as thin film and heterointerface, the 
parameter-free description of correlated systems has be¬ 
come more imperative. In this context, QSGIT study 
of these classical Ad systems can be of significant impor¬ 
tance, and its ability and limitation should be carefully 
investigated. 

Our calculations show that the effect of electronic cor¬ 
relations beyond LDA and GGA is reasonably well cap¬ 
tured by QSGkF procedure, and the noticeable band¬ 
width reduction is observed. For the paramagnetic 
SR0214, the effect of GW self-energy is clearly distinc¬ 
tive from the other many-body calculation {e.g., DMFT) 
in that the spectral weight transfer and the LHB feature 
do not appear in the QSGVF result. For the ferromag¬ 
netic SROII3, QSGkF band structure is quite similar 
with the results of DFT-I-I/ and the half-metallic band 
structure is reproduced. 


II. COMPUTATION METHOD 

A. Quasiparticle self-consistent GW method 

In QSGVF, the self-energy S(r, r',w) is calculated 
within the GW approximation, and Hq (the noninteract- 



FIG. 1: The crystal structure of (a) SR0214 and (b) SR0113. 
The green, gray (octahedra), and red spheres represent Sr, Ru 
(RuOe cage), and O atoms, respectively. 


ing Hamiltonian describing quasiparticles or band struc¬ 
tures) and W (dynamically-screened Coulomb interac¬ 
tions between quasiparticles within random phase ap¬ 
proximation (RPA)) is updated self consistently |55H55] . 
The static nonlocal one-particle exchange-correlation po¬ 
tential U’“’(r,r') is generated from E(r,r',w) as 


1 

= - dwRe[E{uj)]S{uj - H°) + c.c. 

^ J — OO 

I / \ / ; I R-6[^(^i) + ^(^t)] I ; \ / ; 


where and refer to the eigenvalues and eigenfunc¬ 
tions of iLo, respectively, and Re[E(e)] is the Hermitian 
part of the self-energy [SS] [55]. With this one can 
solve a new static one-body Hamiltonian iLo, and con¬ 
tinue to apply GW approximation until self-consistency 
is achieved. The distinctive feature of QSGVF compared 
with DFT-bC/ and DMFT is that it does not require any 
ad hoc parameters. Previous studies, ranging from atoms 
[59] . semiconductors [Ml Ea [701 El] to the various 3d 
transition metal oxides [MIIS3EH71] and 4/ electron 
systems m have demonstrated its capacity to investi¬ 
gate many different types of correlated materials. 


B. Computation details 

Our implementation of QSGIT in ‘ecalj’ code E5] is 
based on the ‘augmented plane wave (APW) -I- muffin- 
tin orbital (MTO)’ method, designated by ‘PMT’ [77] . 
The accuracy of this scheme was proven to be satisfactory 
A key feature of this scheme for QSGVF is that the 
expansion of can be made with MTOs, not APWs, 
which enables us to make the real space representation of 
at any k point. In contrast to the previous approach 
based on FP-LMTO (full potential linearized muffin-tin 
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FIG. 2: The calculated (a) total DOS and (b-d) PDOS for 
SR0214 by LDA (black dotted lines) and QSGIF (red solid 
lines). The spectral function A{uj) calculated by QSGIT self¬ 
energy is plotted in (e) (green solid line) along with the non- 
interpolated total DOS (red solid line). The Fermi level is set 
to zero. 


orbital) [67| . our scheme is free from the fine tuning of 
MTO parameters. 

We used the in-plane and out-of-plane lattice constant 
of 3.87 A and 12.73 A, respectively, for the body-centered 
tetragonal SR0214. For SR0113, the cubic perovskite 
structure is considered as shown in Fig. [^b). The pseu- 
docubic lattice constant of the orthorhombic structure is 
3.93 A [7H]. We used 6x6x6 [73] and 8x8x8 k points 
for the self-energy calculation in the first Brillouin zone 
of SR0214 and SR0113, respectively. For cubic SR0113 
we also performed DFT-l-17 calculations using OpenMX 
code IHO] to make a comparison with previous studies. 
Both LDA [81] and GGA [82] functionals were used in 
combination with Dudarev form of DFT-I-G [83] [84]. 
16x16x16 k points were adopted for DFT-I-C7 calcula¬ 
tions. 



FIG. 3: The calculated band dispersion and the Ru-d PDOS 
for SR0214 by (a) LDA and (b) QSGW. The red, cyan, blue, 
green and gray lines refer to the Ru dyz, dzx, dz. 2 _y 2 
and d ^2 states, respectively. The yellow lines represent the 
O and Sr states (not shown in DOS). The thickness of the 
bands corresponds to the amount of the corresponding orbital 
character. The red and blue arrows represent the bandwidth 
of each band. The Fermi level is set to zero. The insets 
highlight the band dispersion at around —0.6 eV along F to 
M (7r,0,0) by (a) LDA and (b) QSGVF. A sudden change of its 
slope at (0.557r,0,0) is observed only in QSGIF (red arrow). 
This feature is reflected in PDOS (see Fig. [^b)). 


III. RESULT AND DISCUSSION 

A. Electronic structure of Sr 2 Ru 04 


Our LDA result is in good agreement with the previous 
calculations [21H33] . see Fig. (black dotted lines) and 
Fig. The van Hove singularity (vHS) is located at ~ 60 
meV above the Ep and the density of states (DOS) at 
Ep is N{Ep) ~ 3.34 states/(eV f.u.). The Ru-t 2 g state 
is dominant near Ep and is hybridized with 0-2p. From 
Fig. ga), (c), and (d), one can see that the total DOS 

in the range of —4-2 eV is largely attributed to the 

apical 0-2p. The calculated band dispersion is presented 
in Fig. I^a). The bandwidth of the 2D-like d^y state is 
about two times larger than that of dy^px states. The 
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FIG. 4: (a-c) The PDOS for SR0214 calculated by LDA 
(black dashed lines), LDA+DMFT [3^ (blue dotted lines) 
and QSGIF (red solid lines). The Cg states in the previ¬ 
ous LDA and DMFT study were identical and are not 
shown here, (d) The integrated d DOS defined as N{E) = 
n{e)de, where n(e) is the DOS at energy e. Each value 
was normalized to have an equal number of states at E = oo. 
Ep is set to zero. 


nonbonding 0-2p state is located in —4-2 eV as can 

also be seen in Fig. [^d). 

The QSGVF DOS is presented in Fig. (red solid lines) 
|86| . While the overall shape of the DOS is not much 
different from the LDA results, some differences are ob¬ 
served. The two vHS located at ~ -1-60 and ^ -1-260 
meV in LDA [24U2^ become closer in QSGLF shifting to 
~ -1-70 and ^ -1-220 meV, respectively (inset in Fig.j^b)). 
Both the bonding (at ^ —6 eV) and the antibonding 
(around the Ep) part of Ru-4d states becomes flatter in 
their dispersion compared to the LDA result. The t 2 g 
bandwidth is slightly reduced (see Fig. and|^. An 
interesting feature is found in the dxy band dispersion 
along F to M. The different dispersion at ^ —0.6 eV (see 
Fig. i insets) is related to the more pronounced peak 
in the QSGIF PDOS (Fig. j^b)). The 0-2p states are 
shifted to the lower energy regions in QSGFF (Fig. [^c)- 
(d)). The Sr-4d state is further pushed up to ^ 4 eV (not 
shown). 

The further details of Ru-t 2 g state is of particular in¬ 
terest as the LHB feature at around —3 eV has been a 
subject of debate [IHl 031 H5H35] . Fig. compares the 
calculated PDOS by LDA, LDA-bDMFT, and QSGVF. 
Note that the states in the range of —4 -1 eV is no¬ 

ticeably larger in LDA-I-DMFT than LDA and QSGFF 
especially for dy^^zx- It is related to the spectral weight 
transfer from the Tveai-Ep regions to the LHB. This is 
one of the main features of DMFT calculations and was 
the main point of the previous debate between LDA, 
DMFT [ 4511471148 ] and the XPS experiment m, while 
The QSGFF result in this region is similar with the LDA 


rather than the DMFT. Namely, QSGVF procedure does 
not capture the dynamic ‘Mott-Hubbard physics’ well 
as in DMFT while it still takes some correlation effect 
into account as reflected in the bandwidth reduction 
and the mass enhancement. This can also be seen in 
Fig. [^a). One can notice the shoulder-like LHB feature 
in XPS spectrum at ~ — 3 eV being consistent with 
LDA-I-DMFT, while both LDA and QSGVF give the less 
states in this energy region. The theoretical spectra in 
Fig. [^a) were broadened as in Ref. [45]. Not surpris¬ 
ingly, this LHB feature of DMFT is more pronounced in 
dyz,zx states than in dxy because of the narrower band¬ 
width. This point is highlighted in the integrated DOS 
(Fig. I^d)). While the total number of Ru states (the 
integrated value up to Ep) is basically the same in all 
three calculations, the DMFT value gets increased in a 
wider energy range. On the other hand, LDA and QSGVF 
values are much more rapidly increased in the narrower 
energy range of —0.5 < e < 0 eV. Since the QSGVF band¬ 
width is narrower than LDA, this increase is more pro¬ 
nounced in QSGVF. Our result demonstrates the charac¬ 
teristic feature of QSGVF to take into account of electron 
correlations distinctive from DMFT especially regarding 
the incoherent states at —3 eV. 

The Ru-eg bands are also affected by QSGVF self¬ 
energy (see Fig. [^. The slight up-shift of the Cg anti¬ 
bonding bands is consistent with the previous study of 
SrVOa |S7|. As a result, the dx 2 -y 2 band does not touch 
the Ep even at F and Z point (see Fig. |^. The 0-2p 

nonbonding state is located at —3.5-2.5 eV in QSGVF 

(Fig. i, which is in better agreement with the recent 
ARPES result [55] than LDA. 

The overall shape of the calculated Fermi surface by 
QSGVF (see Fig. ^ is consistent with the LDA 
dHvA I^ISlj and ARPES [37VI33] . The inclusion of spin- 
orbit coupling as a perturbative correction within QSGVF 
induces relatively small modifications at the band cross¬ 
ing points (not shown) jSSj. The calculated effective 
mass at each k direction is shown in Fig. [^ with color 
plot where m*/TOLDA is estimated simply by taking the 
derivative of band dispersions. Along the F to X ( 7 r, 7 r, 0 ) 
line, the mass enhancement by QSGVF is about 15% 
(TO*g/TOLDA ^ 1-15) for dxy and 22% (m*^ ,,,,,/tolda ^ 
1.22) for dyz^zx, respectively. For the F to M ( 7 r, 0 , 0 ) line, 
m*y/mLDA — 0.83 and — 0.78. 

It is instructive to compare the QSGVF result with the 
DMFT calculations and experiments while the DMFT 
effective mass strongly depends on the U and J values. 
The value reported by Mravlje et al. with U = 1.7 
and J — O.Q - 0.1 eV is m* jmppix — 1-7 for both dxy 
and dyz^zx [!]• This value is signihcantly smaller than 
the other DMFT calculations, e.g., by Pchelkina et al., 
m*ImpQA — 2.6 and 2.3 for dxy and dyz^zx, respectively 
(with U = 3.1 and J = 0.7 eV obtained from con¬ 
strained LDA) |3^, and by the same group but with 
U = 2.3 and J = 0.2 eV, to^/tulda is — 2.3 and 2.0 
for dxy and dyz,zx, respectively |2|. These DMFT val¬ 
ues were calculated from mDMFx/wLDA = ^dmft = 
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FIG. 5: The experimental and calculated spectra of (a) SR0214 and (b) Ru-4d states of SR0113. In (a), XPS spectrum 
(yellow dots) by Yokoya et al. [40| are plotted along with Gaussian-broadened DOS by LDA (black dashed line), LDA+DMFT 
(blue dotted line) [45], and QSGIF (red solid line). In (b), PES spectrum (yellow dots) by Takizawa et al. |64j obtained from 
orthorhombic SR0113 at room temperature are plotted with Gaussian-broadened DOS of cubic SR0113 (LDA and QSGVF) 
and orthorhombic SR0113 (LDA-I-DMFT [flip. Only Ru-4d states are plotted in (b). The intensities were normalized in the 
area under curves. 



FIG. 6: The Fermi surface of SR0214 at = 0 plane. 
The color represents the calculated effective mass defined as 
W,QSGw/mLDA. 


[1 — 9t.jReSDMFT(w)]w=o- Our estimation from QSGFF 
self-energy yields ^qsgw — 

dyz,zx [Hni which is not much different from the DMFT 
values of 1.7 calculated by Mravlje et al. with U = 1.7 
and J = 0.0 - 0.1 eV |2|. Note that U = 2.3 and J = 0.4 
eV were suggested to be reasonable |5] in comparison 
to dHvA measurements 13011231, and that the experi¬ 
mental values are larger than the QSGFF. The early 
ARPES reports to*/tolda — 2.5 m and the more re¬ 
cent measurement by Iwasawa et al. is m*j^/mLDA — 
3.7 and m *2 2 a;/mLDA — 2.0 |55|. There are two re¬ 
ports from the dHvA experiment; m*y/mLDA — 4.1 
and TO* 2 _j;a;/wLDA — 3.3 by Mackenzie et al. [30], and 
mlylmi.DA — 5.5 and — 3.4 by Berge- 

mann et al. |32j . 


B. Electronic structure of cubic SrRuOa 

The LDA result for cubic SR0113 is in good agree¬ 
ment with the previous studies [TT] [S0H23] (see Fig. 


and Fig. [^. The calculated magnetic moment of ^ = 
1.28^B/f-u. also reasonably well agrees with the litera¬ 
ture values of /r = l.OO/xe/f.u. (calculated by VASP) and 
fjL = 1.26/iB/f.u. (calculated by SIESTA) [T7]. The Sr-dd 
state is located above ^ 4 eV (not shown) and the non¬ 
bonding state of 0-2p is in —4- 2 eV. The antibonding 

Ru-t 2 g character dominates the neai-Ep region and the 
bonding complex is located at —8-4 eV. A clear split¬ 

ting between the up and down spin DOS is noticed and is 
responsible for the ferromagnetism in this material. The 
effect of orthorhombic distortion (not taken into account 
in our calculation) has been investigated previously. For 
example, Rondinelli et al. m showed that the result 
from orthorhombic structure is similar to the cubic case 
other than the slightly reduced exchange splitting at F 
point and bandwidth reduction by ~ 0.35 eV for t 2 g, 1.5 
eV for eg, and 0.6 eV for O 2p. 

Several distinctive features are found in the QSGIE re¬ 
sults. First, the notable bandwidth reduction is observed 
(see Fig. and Fig. [^. The majority spin bandwidth is 
reduced by ^ 0.7 eV (from ~ 3.2 eV (LDA) to ^ 2.5 eV 
(QSGIF)) and the minority spin bandwidth is by ~ 0.6 
eV (from ~ 3.2 eV (LDA) to ~ 2.6 eV (QSGVF)). The 
exchange splitting is enhanced to be ~ 1.2 eV which is 
significantly larger than the LDA value of ~ 0.5 eV. We 
found the naive comparison of the effective mass for LDA 
and QSGFF based on the k-derivative can be misleading 
in SR0113 due to the different Fermi wave vectors caused 
by the enhanced exchange splitting in QSGIF. The cal¬ 
culated ^QgGvv down spin is 1.33 and 1.68, 

respectively. The effective mass of toqsgw/wlda = 1-26 
estimated from the bandwidth ratio is reasonably well 
compared with some of DMFT results while the DMFT 
values show significant deviations ranging from 1.1 to 4.5 
(59] IMHM] . The specific heat measurements report the 
mass enhancement of TO*/mLDA = 3.7 [SD] and 4.5 [61j . 
The discrepancy reflects not only the limitation of RPA- 
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FIG. 9: The calculated Ru t 2 g and Cg PDOS of cubic SR0113 
by LDA+f/ {UeB = 2.5 eV) (blue dashed lines), GGA+17 
{Ues = 2.0 eV) (green dotted lines) and QSGIT (red solid 
lines). The Fermi level is set to zero. The positive and nega¬ 
tive DOS represent the up- and down-spin parts, respectively. 


FIG. 7: The calculated (a) total DOS and (b-c) PDOS of 
cubic SR0113 by LDA (black dotted lines) and QSGVF (red 
solid lines). The Fermi level is set to zero. The positive 
and negative DOS represent the up- and down-spin parts, 
respectively. 



FIG. 8: The calculated band dispersion and the Ru d-orbital 
PDOS for cubic SR0113 by (a) LDA majority spin band, (b) 
LDA minority spin band, (c) QSGVF majority spin band, and 
(d) QSGVF minority spin band. The red, cyan, blue, and gray 
lines refer the Ru d^y, dyz, dzx, and Cg states, respectively. 
The yellow lines represent the O and Sr states (not shown in 
DOS). The thickness of the bands corresponds to the amount 
of the corresponding orbital character. The Fermi level is set 
to zero. 


QSGLF correlations but also the effect of orthorhombic 
distortion isiMiiini- 

Together with the bandwidth reduction and the en¬ 
hanced exchange splitting, the QSGIT electronic struc¬ 
ture becomes half-metallic with a gap in the majority 
spin state. Note that this QSGLF band structure is 
quite similar with DFT-I-C/ result nuns] as shown in 
Fig.§ The calculated magnetic moment of QSGFF is 
/i = 2.0y(B/f.u. which is significantly larger than the 
LDA value and comparable with the DFT-I-C/ results 
[IMZj. While we only considered the cubic structure, 
the effect of orthorhombic distortion is basically to en¬ 
hance the on-site correlations as reported by the previ¬ 
ous study [6^ [90] . The experimental verification of the 
half-metallicity may not be easy because of the large mag¬ 
netic fields required to overcome the magnetic anisotropy 
m- This possibility in SR0113 has been actively dis¬ 
cussed based on DFT-I-C/ calculations [T5lll6| and hybrid 
functional [SSj. Therefore our QSGLF result adds a new 
promising aspect toward this direction. 

The detailed comparison with DFT-I-C/ result is given 
in Fig. The optimized values of Ues are favorably 
compared with a recent constrained RPA result of C/eff = 
2.1 eV although it is calculated from the orthorhombic 
structure [T^. For comparison one can recall that SIC- 
LDA result is reproduced by C/eff = 1 eV in the case 
of orthorhombic structure of SR0113 m- As for the 
DMFT calculations, several different choices of U and J 
have been made as in the case of SR0214. In terms of 
Uee=U — J, it ranges from 1.75 to 2.4 eV 

Several photoemission spectroscopy (PES) experi¬ 
ments report the LHB-like feature at ~ —2 eV [60l 
iSIllSl]. The PES spectrum of orthorhombic phase |64] 
are plotted along with the calculated DOS by LDA, 




















































































7 


LDA+DMFT [51], and QSGIF in Fig.[5];b). The LHB- 
like peak near — 1.5 eV is observed in the PES spectra, 
while not in LDA and QSGVF. It has been argued that 
this LHB-like feature is not related to the magnetic fluc¬ 
tuations or the orthorhombic structural distortion [M] , 
LDA-I-DMFT calculation of orthorhombic SR0113 with 

U = 3.5 and J = 1.75 eV predicts the LHB at —3-1 

eV |5T] (Fig.j^b)). As in SR0214, LHB is the character¬ 
istic feature of DMFT distinctive from QSGVF (also from 
DFT-I-G) while some DMFT results do not seem to sup¬ 
port this feature [S5]. The absence of this state in QSGkF 
and DFT-I-I/ (see t 2 g states of Fig. 1^ indicates that it 
is related to the dynamic aspects of Ru-4c? correlation 
hybridized with 0-2p as in the case of Zhang-Rice band 
in cuprates. Also, regarding the half-metallic band struc¬ 
ture, DMFT calculations do not seem to give a consistent 
prediction |9TJ |93| . These issues that is likely related to 
the intrinsic parameter dependency require further inves¬ 
tigations. 

IV. SUMMARY 

Using QSGkF calculations, we re-investigated the elec¬ 
tronic structure of SR0214 and SR0113. Without any 
ad hoc parameter or the ambiguity related to the double¬ 
counting and downfolding issues, some of the important 
features of electron correlations were reasonably well cap¬ 
tured such as the bandwidth renormalization and the ex¬ 
change splitting. In the case of SR0113, QSGkF result is 


in good agreement with DFT-l-C/ in a reasonable range of 
U and J parameter. While the QSGIU shows the limita¬ 
tion in describing the detailed features such as the Ru-dd 
spectral weight transfer to LHB, it can be improved in 
combination with other techniques as reported in recent 
studies mM [95] . Our result sheds new light on the 
possibility and the limitation of the first-principles elec¬ 
tronic structure calculations of the moderately correlated 
transition-metal oxides systems. 
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